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(10 x 10) SWNT consisting of 400 atoms with 20 layers is simulated under tensile loading using 
our developed O(N) parallel TBMD algorithms. It is observed that the simulated carbon nanotube 
is able to carry the strain up to 122% of relaxed tube length in elongation and up to 93% for 
compression. The Young's modulus, tensile strength and Poisson ratio are calculated and the values 
found are 0.311 TPa, 4.92 GPa and 0.287, respectively. Stress-strain curve is obtained. The elastic 
limit is observed at the strain rate of 0.09 while the breaking point is at 0.23. The frequency of 
vibration for the pristine (10 x 10) carbon nanotube in radial direction is 4.71 x 10 3 GHz and it is 
sensitive to the strain rate. 
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Single walled carbon nanotubes (SWNT) are observed 
to have high flexibility, strength and stiffness, very sim- 
ilar to those of individual graphene sheets. Eventhough 
direct measurements of such mechanical properties are 
difficult to perform due to the nanosizes involved; along 
with the developments in instrumentation, production, 
processing and manipulation techniques, the measure- 
ments of the elastic moduli of carbon nanotubes has be- 
come possible. On the other hand, this extreme small 
size is very suitable for performing atomistic simulations. 
Both the experimental^-^ and theoretica l 3 15 studies 
have shown that SWNTs, and SWNT ropes are promis- 
ing low- weight, high-strength fibers for use as reinforcing 
element in composite materials. It is also predicted that 
SWNT's can sustain large strains in the axial direction 15 . 
The axial Young's modulus values range from 200 GPa 
to 5.5 TPa in literature. 

The carbon nanotubes are known to have honeycomb 
structure^ and three C-C bonds formed between a car- 
bon atom and three adjacent carbon atoms in the unit 
cell can be classified as two kinds according to the rela- 
tions of their spatial orientations; one is perpendicular to 
tubular axis and the other is not perpendicular to tubu- 
lar axis in armchair tube. In this study, an armchair 
(10 x 10) carbon nanotube consisting of 400 atoms with 
20 layers is simulated under tensile loading using our de- 
veloped O(N) parallel TBMD programii^. Two steps 
are followed; firstly, the tube is annealed at simulation 
temperature for 3000 MD steps (each time step used in 
the simulation being 1 fs). The variation of total energy 
and some physical properties such as radial distribution 
function, atomic coordination number, bond-angle dis- 
tribution function, and bond-length distribution func- 
tion are given. The second step is to apply strain on 
the tube. Several groups have proposed different proce- 
dures; such as shifting the end atoms along the axis (i.e. 
z-direction) by small steps^, reducing radial dimension 
while the nanotube is axially elongatedAS, pulling in the 



axial direction with a prescribed strain rate and follow- 
ing each step of pulling by some additional MD steps in 
order to relax the distorted structure^ and finally us- 
ing hydrostatic pressure exerted on the walls of SWNT 
by the means of encapsulating H2 molecules inside the 
tube and the wall of the tubeii In our study, the tensile 
strain is applied by the reduction or enlargement of the 
radial dimension while the nanotube is axially elongated 
or contracted. Throughout this procedure volume of the 
tube is kept constant. Zhou et aliiS, has investigated the 
mechanical properties of SWNT with the same procedure 
using a first-principle cluster method within the frame- 
work of local density approximation. We further simu- 
lated the deformed tube structure (under uniaxial strain) 
for another 2000 MD steps (time step is chosen again as 
1 fs) to understand the strain mechanism. Strain is ob- 
tained from e = ( L ^ ^ , where Lq and L are the tube 

lengths before and after the strain, respectively. Several 
strain values are applied to a pristine tube to study the 
strain rate. Simulations are performed at room temper- 
ature and periodic boundary condition is applied along 
the axial direction. 

The binding energy curve of a carbon nanotube as a 
function of the strain along the tubular axis is given in 
Figure^ We observe an asymmetric pattern for the cases 
of elongation and compression. The tube does not have 
a high strength for the compression as much as for elon- 
gation. This might be due to the dominant behavior 
of repulsive forces in the system under uniaxial strain. 
This figure indicates that the remarkable elastic proper- 
ties under large strains are caused by nonparabolic strain 
energy. We observed that the carbon nanotube is able to 
carry the strain up to 122% of the pristine tube length 
in elongation and up to 93% of pristine tube length in 
compression. 

The variation of total energy of the deformed system 
during MD simulation for the strains of 0.22 in elongation 
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and 0.07 in compression are given in the Figure|2J In the 
graphs, first 3000 MD steps is for the equilibration of 
the carbon nanotube and the next 2000 MD steps shows 
the variation of the total energy of the tube during the 
simulation under the applied uniaxial strain. It is seen 
that the tube under these strain rates is able to sustain 
its structural stability.For high strain rates the changes in 
the radial distribution function, bond-length distribution 
and bond-angle distribution are given in Figs. I2HSI It is 
seen in the graphs that the two-third of the bond lengths 
and one-third of the bond angles increases (decreases) 
when tube is elongated (contracted), as expected. Bond 
angles and bond lengths are the two important factors 
that control the deformation. The effect of the strain on 
the bonds is such that it alters the angles between two 
neighboring carbon bonds and changes the lengths of the 
C-C bonds. 

Increasing the strain beyond these points results in dis- 
integration of atoms from the carbon nanotube. The ge- 
ometrical structures, and the behavior of the total energy 
for the strains 0.23 (elongation) and 0.08 (compression) 
show that the elongated tube dissociates by starting from 
the middle like a zipper while the compressed tube starts 
to dissociate from the ends of the tube. Each peak in 
Figure El represent disintegrations of the atoms from the 
tube. 

The elastic constants are calculated from the second 
derivative of the energy density with respect to various 
strains. To obtain the stress-strain curve, the cross- 
section upon which the resulting forces act is needed to 
be estimated. The cross-sectional area of a nanotube is 
ambiguous in definition 16 . If a circular cylindrical shell 
is considered around the surface of the nanotube, then 
the surface area of the cross-section s, is defined by 
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where R stands for the radius of the SWNT and 5R for 
the wall-thickness. It should be noted that different wall- 
thickness values were used by several groups 5 - 8 . In Ref. 7 , 
<5i?=3.4 A (measured interwall distance in the Multi Wall 
Nano Tube) was used, while in Ref. 8 5R=1.7 A (taken 
as the van der Waals radius for Carbon) was accepted , 
it is also accepted as 0.66 A (in the tt orbital extension) 
in Ref£ and as the whole cross-sectional area of the tube 
in the RefA We defined the thickness of SWNT shell 
as SR = 3.4 A. The stress-strain curve obtained from 
this study is given in the Figure [7] It is seen in the figure 
that the elastic limit is at the strain value of 0.09. Beyond 
the elastic limit, the stress-strain curve departs from a 
straight line. Hence, its shape is permanently changed. 
The breaking point is observed at the strain rate of 0.23. 
The Young's modulus is determined as the slope of the 
stress-strain curve. Our calculated value of the Young's 
modulus of the (10 x 10) Carbon nanotube is 0.311 TPa. 

Theoretical tensile strength is defined as the maximum 
stress which may be applied to the material without per- 



turbing its stability. It can be given as 
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where e is the stress, £j is the maximum stress in the 
system, and s is the surface area of cross-section. Our 
calculated value is 4.92 GPa, which is larger than that 
of Carbon fibers (2.6 GPa)i^, but less than the in-plane 
tensile strength of graphite (20 GPa)Si. 

Another mechanical property of interest is the Poisson 
ratio, defined by 
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where R is the radius of the tube at the strain e, and R eq 
is the equilibrium (zero strain) tube radius. The Poisson 
ratio measures how much the tube contracts (expands) 
radially when subject to a positive (negative) axis strain 
e. The corresponding value found in this study is 0.287. 

Another interesting phenomenon we observed in sim- 
ulation is the vibration of SWNT in radial direction. In 
Figure |S1 the average radius of the pristine (10 x 10) 
SWNT as a function of MD steps is given. The frequency 
of vibration can be evaluated from the figure and has the 
value of 4.71 x 10 3 GHz. The variations of radius for the 
strains 0.22 and -0.07 are also given in the Figure [5] It 
is found that the strain is effective on the vibration fre- 
quency. Increasing strain on the tube structure results in 
the decrease for the frequency of vibration. It is ranged 
from 2.94 x 10 3 GHz to 4.41 x 10 3 GHz with decrease 
by the increasing strain rate and has a mean value of 
3.70 x 10 3 GHz. 

We have thus determined the elastic properties of a 
(10 x 10) carbon nanotube under tensile loading and 
found the Young's modulus, tensile strength, Poisson ra- 
tio and frequency of vibration to have the values 0.311 
TPa, 4.92 GPa, 0.287 and 4.71 x 10 3 GHz, respectively. 
Several groups have reported a wide range of values for 
the corresponding properties by using various theoretical 
and experimental techniques. 

The Young's modulus values given by different re- 
searchers range from 0.200 TPa to 5.5 TPa 14 . The fol- 
lowing reasons may be given for this variety of results: 

• The different values are used for the wall- 
thickness^. 

• Different procedures are applied to represent the 

• The curvature effect of nanotubes was neglected^, 
or not 13 . In Ref. 7 , it is concluded that the elastic 
moduli of nanotubes, (SWNT and MWNT) were 
insensitive to geometrical structure while it is sug- 
gested in Refji 3 - that Young's modulus slightly de- 
pend on the tube diameter. On the other hand, the 
variation of Young's modulus as a function of tube 
radius is also reported^ 6 -. 
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• Accuracy of methods: first-principle methods are 
more reliable^jiSji 2 "^ with comparison to the em- 
pirical potentials 3 *^. 

• The strain rates that Young's modulus was calcu- 
lated are either differenfcii or not pointed oul* 7 ^. 

• Difference in the tube lengths: although the pe- 
riodic boundary condition is applied for the most 
cases, finite size effects might be still important. 

Our result is in the range mentioned above. It emphasizes 
the high Young's modulus and high strengths of carbon 
nanotube. The strain at tensile failure for SWNTs was 
predicted to be as high as 0.4C& This is a tensile strength 
of 400 GPa would be expected for SWNTs if one used the 
in-plane Young's modulus of graphite, ~ 1 TPaii. How- 
ever, such a high tensile strength has not been justified 
by experiments. In this study, it is found that the elastic 
limit is at the strain rate 0.09 and beyond this point tube 
becomes permanently changed. In Refi-a, it is reported 
that for strain values greater than 0.10, the tube becomes 
softened. They also estimated the strain at failure for the 
SWNT as 0.17 whereas it is found as 0.23 in this study. 
The procedure for describing the strain in their work is 
to apply hydrostatic stress to the tube wall. On the other 
hand, Tight-Binding electronic calculations reported by 
Ozaki et al& revealed a strain as high as 0.30. 

The calculated and measured tensile strength also 
varies in value. In Ref^, it is reported as 6.249 GPa by 
a result of a first-principle study while it has the value 
of 62.9 GPa for the perfect 5x5 SWNT under hydrostatic 



pressure!. 4 . Another MD simulation by using the multi- 
body potential function of embedded atom method re- 
ports the tensile strength as 9.6 GPa£. On the other 
hand, it is reported as 3.6 GPai and as ranged from 13 
to 52 GPa— in the experimental studies. The value found 
in this study is 4.92 GPa and seems to comparable with 
the experimental and theoretical results. 

The calculated Poisson ratio is 0.287 and in good agree- 
ment with the available reported values which are 0.278 7 , 
and 0.32 10 . The evaluated frequency of vibration for the 
pristine (10 x 10) Carbon nanotube is 4.71 x 10 3 GHz 
which is very close to the value obtained from the exper- 
iment 4.94 x 10 3 GHaSi and almost same with the value 
reported in the MD simulation by using a bond-order 
potential 12 . In Ref. 12 , it is reported that the frequency 
of vibration is insensitive to the strain rate and the fre- 
quency of vibration is identified as self-vibration. We 
have found that it is not constant and increasing the 
strain rate decreases the vibration frequency. 
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FIG. 1: Total energy curve as a function of strain e 
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